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Abstract 

We calculate the exact eigenvalues of the adjoint scalar fields in the massive vacua of 
J\f = 1* SUSY Yang-Mills with gauge group SU(N). This provides a field theory prediction 
for the distribution of D3 brane charge in the AdS dual. We verify the proposal of Polchinski 
and Strassler that the D3-brane's lie on a fuzzy sphere in the supergravity limit and determine 
O ! the corrections to this distribution due to worldsheet and quantum effects. The calculation 
also provides several new results concerning the equilibrium configurations of the iV-body 
Calogero-Moser Hamiltonian. 
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1 Introduction 



The N = 1* deformation of Af = 4 SUSY Yang-Mills theory with gauge group SU(N) 
provides a fascinating example of a non-conformal gauge theory where a direct comparison 
of the AdS dual with field theory results is possible. Specifically, as this is a theory with 
M = 1 SUSY, there are special observables which depend holomorphically on the parameters 
and can be computed exactly for all values of the marginal coupling r = Airi/g 2 + 9/2% and 
N. In the limit of large 't Hooft coupling, A = g 2 N/47r — > oo, the results can be compared 
with a supergravity computation in the IIB dual proposed by Polchinski and Strassler JT7I . 



The M = 1* theory contains three complex scalar fields <£>i, $2 and $3, in the adjoint 
representation of G = SU(N). After a rescaling, these fields appear on an equal footing 
and may be rotated into each other by an SO (3) subgroup of the 577(4) R-symmetry of the 
M = 4 theory. In the following we will study the eigenvalues of these N x N matrices in 
super symmetric vacuum states of the theory, or equivalently, the vacuum expectation values 
of the N — 1 gauge- invariant operators Uk = Tr^^ k for k = 2,...,N. Here, and in the 
following, $ denotes any one of the three $j and Afe, with k = 0, 1, . . . , N — 1, denotes its 
eigenvalues. The main result of this paper is an exact formula for the large- N distribution 
of these eigenvalues in each massive vacuum state of the theory. As we explain below, the 
desired result may be obtained by considering the equilibrium configurations of a certain 
classical integrable system. In the remainder of this introductory section we will give this 
result and discuss its physical interpretation. Details of the calculations are provided in 
subsequent sections. 



We begin by reviewing the classical vacuum structure of the theory. The classical F-term 
vacuum equations are solved when the three adjoint scalars (appropriately rescaled) obey 
an 577(2) algebra $-,•] = In fact, as first explained in 0, the classical theory has 

one SUSY vacuum for each iV-dimensional representation of 577(2), with the VEVs of the 
$i being given by the ant i- Hermit ian generators acting in this representation. An arbitrary 
representation of dimension N can be decomposed as the sum of irreducible representations 
whose dimensions add up to N. The total number of classical vacua is therefore equal to the 
number of partitions of N. The unique irreducible representation corresponds to a vacuum 
where the gauge group is completely broken by the Higgs mechanism and the spectrum is 
massive. We will refer to this as the 'Higgs vacuum'. In vacua corresponding to reducible 
representations, a non-trivial subgroup of SU(N) is left unbroken and there are massless 
particles at the classical level. 



This correspondence between classical M = 1* vacua and SU(2) representations has an 
interesting realization in weakly coupled IIB string theory where the undeformed M = 4 
theory with gauge group U(N) lives on the worldvolume of N parallel D3 branes. The 
M = 1* deformation corresponds to turning on the flux of the two-form antisymmetric tensor 
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field of the Ramond-Ramond sector of the IIB theory. The background RR flux causes the 
D3 branes to polarize to spherically wrapped D5 branes by the Myers effect ||. As there is 
non-zero D3-brane charge on the two wrapped dimensions of the fivebrane, the corresponding 
worldvolume theory is non-commutative. In field theory, the non-commutivity is realized by 
the 577(2) algebra described above satisfied by the VEVs of the adjoint scalar fields, $j. The 
N- dimensional irreducible representation corresponding to the Higgs vacuum is the familiar 
matrix representation of the fuzzy sphere. Reducible representations correspond to multiple 
concentric fuzzy spheres in an obvious way. 



With the standard normalization of the generators, the eigenvalues of each of the adjoint 
scalars in the Higgs vacuum is uniformly distributed on the imaginary axis between ±i(N — 
l)/2 with unit spacing. Explicitly, the N eigenvalues are = i(N — 2k — l)/2 with k = 
0,1, N — 1. In Section 2.2, we will derive an exact quantum mechanical formula for the 
eigenvalues A& in the Higgs vacuum (see equation (p8|). In the large- N limit we replace 
k/N by a continuous variable < x < 1. The range of the eigenvalues determines the 
radius of the sphere on which the D3 branes lie. Equivalently, it is also determined by the 
expectation value of the lowest non-zero condensate u 2 which has already been studied in 



detail in earlier work [II], [13], [12|]. For the present purposes it is convenient to normalize the 
range of the eigenvalues to unity, defining a normalized eigenvalue distribution \{x). (The 
natural normalization as we later show involves the dual 't Hooft coupling and the modular 
weight.) The large- N limit of the semiclassical result therefore yields a uniform distribution 
X(x) = 1 — x/2. In the following we will see that this uniform distribution also emerges in 
the strong-coupling supergravity regime. 



So far we have only described the classical limit of the M = 1* theory. In vacua with 
an unbroken non-abelian gauge symmetry, the low-energy theory is realized in a strongly- 
coupled confining phase with gluino condensation and a mass gap. The vacuum structure of 



the quantum theory is determined by the exact superpotential derived in 11 



W = m im2 m 3 J2V(X a -X b ) (1) 

a>b 

Here, V(X) is the Weierstrass function and rrii, i = 1,2,3, are the masses of the three 
chiral multiplets which we will set to one from now on. The N auxiliary variables X a , with 
a = 1,2, . . . , N, only have a direct physical interpretation when the theory is compactified 
on a circle down to three dimensions. In this context they can be identified with a complex 
combination of the Wilson and 't Hooft loops of the gauge field around the compact direction. 
These variables naturally take values on a torus whose complex structure is the complexified 
coupling constant r, and the exact superpotential, being an elliptic function, is single-valued 
on this torus. The S-duality transformations of the M = 4 theory correspond to modular 
transformations of the torus. 
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Stationarizing the superpotential with respect to the variables N variables X a , we find a set 
of massive vacua labeled by integers p,q and k, with pq = N and k — 0, . . . , q — 1. The total 
number of these vacua is the sum of the divisors of N. For each divisor q of N, the q quantum 
vacua labelled by k = 0, . . . ,q — 1 reduce in the weak coupling limit to the unique classical 
vacuum where the unbroken gauge symmetry is SU(q). This classical vacuum corresponds 
to the reducible representation of SU (2) which is the direct sum of q copies of the irreducible 
representation of dimension p = N/q. The Higgs vacuum corresponds to q = 1, p = N. For 
q > 1, the low energy theory is strongly coupled and gluino condensation occurs, splitting 
the massless classical vacuum into q massive quantum vacua. The theory also has a large 
number of vacua which remain massless at the quantum level but we will not consider these 
here. 



A remarkable feature of the M = 1* theory is that the set of J2d\Nd vacua described 
above are transformed into each other by the S-duality of the underlying Af = 4 theory. 
In particular, the electromagnetic duality transformation r — > — 1/r interchanges the Higgs 
vacuum in which electric charges condense and a confining vacuum in which magnetic charges 
condense. This provides a concrete realization of the duality between Higgs and confining 
phases conjectured many years ago by 't Hooft. This property is manifest in the exact 
results following from the superpotential ([I]). The action of S-duality can also be understood 
in the IIB realization of the theory described above. As above the Higgs vacuum is realized 
on N parallel D3 branes polarized into a spherically wrapped D5-brane. The D3 branes 
are invariant under the SL(2, Z) duality of IIB string theory, but the D5 brane is not. 
The action of S-duality generates an SL{2, Z) orbit of vacuum configurations described by 
(m, n)-fivebranes wrapped on non- commutative two-spheres each carrying N units of D3- 
brane charge. Here m and n are determined in terms of the the integers p, q and k which 



describe the field theory vacua [17 



Polchinski and Strassler studied the AdS duals of these M = 1* vacua. They proposed a 
dual IIB geometry for each vacuum which is asymptotically AdS$ x S 5 near the boundary, 
but is deformed by the presence of the spherically wrapped (m, n)-fivebranes in the interior. 
The precise domain in which the supergravity approximation is reliable depends explicitly on 
the choice of vacuum. We start from the vacuum which exhibits pure magnetic confinement, 
corresponding to the choice p = 1, q = N and k = 0. This vacuum is described by a single 
spherically wrapped NS five-brane. The SUGRA description is valid when the sphere is 
large in string units and this is true as long as the 'electric' 't Hooft coupling, A = g 2 N/An, 
is much larger than one. Away from this limit corrections from IIB worldsheet instantons 
wrapping the S 2 are unsuppressed. Part of the charm of this model is that these corrections 
can be calculated explicitly in the dual field theory. As explained above the remaining vacua 
are related to the magnetic confinement vacuum by S-duality. Specifically, we perform an 
SL(2,Z) transformation, r — > td = (clt + b)/(cr + d) = 4m/gjj + #d/27t where a, b, c, 
d are integers with ad — be = 1 which transforms the wrapped NS five-brane into an (a, c) 
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five-brane. In this vacuum we define a dual 't Hooft coupling Xd = glb^/Air and the classical 
supergravity approximation is reliable whenever Xd » 1 with g D << 1. Away from large 
Xd, corrections coming from (a, c)-strings wrapping the two-sphere become significant. 



In the regime described above, the supergravity dual can be used to compute various field 
theory quantities such as the tensions of BPS domain walls and the VEVs of chiral operators. 
In [1^] , these were successfully compared with the exact field theory results of [[□], |i~3 |. 



However, the particular quantities tested so far correspond to operators of low .R-charge such 
as the gluino condensate and the operator w 2 defined above. These operators correspond to 
massless scalars present in five-dimensional gauged supergravity. One of the main aims of 
this paper is to extend this test to operators of higher R-charge corresponding to Kaluza- 
Klein modes which probe the full ten-dimensional geometry. Specifically, the eigenvalues A& 
determine the VEVs of the operators Uk, which have charge proportional to k under a U(l) 
subgroup of the SU(4) R-symmetry of the M = 4 theory. The corresponding SUGRA states 
live in the tower of spherical harmonics on the five-sphere of the UV geometry. 

The VEVs of the operators can be read off from the asymptotic behaviour of these 
SUGRA fields in the Polchinski-Strassler solution corresponding to each vacuum. Indeed, 



this was the proceedure followed in [12j to compute the VEV of u 2 . However it is simpler to 
note that, as at weak-coupling, the eigenvalues A& determine the position of the distribution 
of D3 brane charge in spacetime. The correspondece between these two points of view for 
the M = 4 theory is discussed in 

0. By calculating the A^ directly in field theory and then taking the SUGRA limit, 
we will be able to check directly that the 3-branes lie on non-commutative two-spheres in 
spacetime as predicted by Polchinski and Strassler. Our main result is a universal formula 
for the large- N eigenvalue distribution which is valid in all the massive vacua 

I oo / e -2%\ D (l-x) e -2nX D {l-l+x) \ 

K X ) = 2 _X + ^^l + e -2n\ D (l~x) ~ I + e -2nX D (l-l+x) J ^ 

Here, as above Xd = 9d-^/^ 1t * s the appropriate 't Hooft coupling in each vacuum. In 
the limit Ad — > oo only the first term on the right-hand side contributes and we find a 
uniform distribution of eigenvalues centered at the origin x = 0. This is consistent with the 
fuzzy sphere structure described above. In particular, we recall that the eigenvalues of the 
generators of SU(2) in an irreducible representation are uniformly distributed around the 
origin. The corrections correspond to the contributions of wrapped (a, c)-strings around this 
sphere [0, 0. 



Our approach to determining the exact eigenvalues A& is based on the connection between 
super symmetric gauge theories and integrable systems first discovered in 0, p| ||. In partic- 
ular, it is noteworthy that the exact superpotential coincides with the potential energy of the 
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elliptic Calogero-Moser system. The latter is an integrable model describing non-relativistic 
particles moving two-body interaction potential determined by the Weierstrass potential. In 
terms of the positions and momenta of the N particles, X a and P a with a — 1, . . . , N, the 
Hamiltonian is, 

N p2 

H=Ef+EP(Ia-I 6 ) (3) 

a=l Z a>b 

Although the positions and momenta are real in the first instance, it is convenient to promote 
them to complex variables. 



In its original form the connection between integrable systems and supersymmetric gauge 
theory involves theories with M = 2 supersymmetry. The theory relevant to the present case 
is the Af = 2 supersymmetric deformation of the M = 4 theory, also known as the M = 2* 
theory. In particular the complex curve which governs the Coulomb branch of this theory 
coincides with the spectral curve of the Calogero-Moser system which plays a key role in its 
explicit integration. As explained in [11], the fact that the exact superpotential of the M = 1* 



theory coincides with the Calogero-Moser Hamiltonian follows from soft breaking to M = 1 
supersymmetry. In particular, the supersymmetric vacua are in one to one correspondence 
with the equilibrium positions of the (complexified) Calogero-Moser system. The results 
given below extend the existing results in the literature about the equilibrium properties of 
the integrable system. In particular, they provide an elliptic extension of the results given 
in m for the trigonometric Calogero-Sutherland model. 



For the present purposes, certain aspects of the correspondence between the M = 1* system 
and the Calogero-Moser system are particularly useful. The integrability of the Calogero- 
Moser system is equivalent to the existence of N conserved quantities or Hamiltonians, 
Hk with k = 1,2, ...N, of which the first two are the total momentum and energy. As 
discussed in the next section, Hamilton's equations of motion for the system (|3]) have a Lax 
Representation and the N conserved quantities are determined by the eigenvalues of the 
Lax matrix L as = Tr^L k . On the other hand, the identification of the Calogero-Moser 
spectral curve with the complex curve of the M = 2* theory reveals that the Hamiltonians 
iffc should be identified with the VEVs of the operators^] = Tr^^ k . This gives us a 
very convenient way to calculate the eigenvalues in each SUSY vacuum: they can be 
identified with the eigenvalues of the Lax matrix in the corresponding equilibrium position 
of the Calogero-Moser system. 



The identification described above is complicated by the problem of operator mixing as 
discussed in [12]. In particular, once masses are turned on, the operator Uj can mix with 

1 The vanishing of u\ in the SU(N) theory corresponds to choosing the center of momentum frame 
Hi = J2a=i Pa = in the integrable model. 
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the other operators Uk with k < j which have lower dimensions. This is particularly clear in 
the M = 2* theory, where the ambiguity is simply that of choosing holomorphic coordinates 
on the Coulomb branch. Similarly, in the integrable system, it is possible to find many 
inequivalent Lax matrices whose eigenvalues correspond to different linear combinations of 
the conserved quantities. Hence there is no unique way to the operators Uk in terms of 
the conserved Hamiltonians. Fortunately this amibiguity turns out to be very mild in the 
large- N limit. As explained in \ Y1\ the ambiguity is vacuum independent and can be fixed by 
determining the mixing coefficients in any given vacuum. These are holomorphic functions of 
the complexified coupling which necessarily have a weak coupling expansion in integer powers 
of the Yang- Mills instanton factor qyu = exp 2nir. These terms vanish exponentially in the 
't Hooft limit, and the mixing coefficients are pure numbers which can be determined by 
a semiclassical calculation. The upshot is that once we have determined the semiclassical 
relation between the eigenvalues of the Lax matrix and the field theory operators in the 
Higgs vacuum we can make definite statements about the large- N distribution of D3 branes 
in the other vacua. 



The paper is organised as follows. In Section 2 we introduce the integrable system and 
determine the eigenvalues of the Lax matrix in the equilibrium configurations. In Section 3, 
we check the prediction of Polchinski and Strassler that the eigenvalues lie on a fuzzy sphere 
in each vacuum. Section 4 we perform a detailed analysis of the results in the confining 
vacuum and in Section 5 we compare with earlier results for the condensate u-i- Another 
approach to calculating the field theory VEVs, based on the complex curve of the Af = 2* 
theory is described in an appendix. 



2 The Calogero-Moser integrable system 

The CM system |7j consists of N non-relativistic particles with position X a which interact 
via the two-body potential 

V = J2'P(X a -X b ) a,b = l...N (4) 

a>b 

where V is the Weierstrass function defined on a torus with periods 2uj\ = 2m and 2uj 2 = 
2mr. The coordinates and momenta are considered complexified and the positions of parti- 
cles are centered so that 

E^ = o (5) 

(1=1 

The system has a non-trivial equilibrium position at, 

X a = i '- a = 1 ... N 6 

N N K J 
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where the N particles are equally spaced along the one-cycle of the torus parallel to the 
period 2u\. As explained in Jll|, this configuration corresponds to the Higgs vacuum of 
the M = 1* theory. Acting by modular transformations, we find a family of equilibrium 
configurations equal in number to the sum of the divisors of N less one where the particles 
are equally space along different cycles of the torus. These correspond to field theory vacua 
in the confining phase. 

The Calogero-Moser system can be described in terms of a pair of Lax matrices, and is 
characterized by N integrals of motion and corresponding angular variables. The Af = 1* 
vacua are described in terms of the N order parameters in the moduli space = Tr$ fc , 
which correspond to the N integrals of motion of the integrable system modulo the potential 
ambiguities described above. The integrals of motions or Hamiltonians for the system are 
the symmetric polynomials of the eigenvalues of the Lax matrix L. 



Hamilton's equations for system of N particles are equivalent to the Lax matrix equation 

1 

iL = [M, L] (7) 



where 



Lab = pJab + i(l - S ab )x(X a - X b ) 
M ab = 5 ab J2z(X a -X c )-(l-5 ab )y(X a -X b ) 



with 



V{u) = x{uY + C 

y(u) = - X '(u) 
. , x"(u) , . 

2( "» = 2^y < 9 > 

where C is a constant. In the present case where the potential V is the Weierstrass function, 
x(u) is given in terms of Jacobi elliptic functions 0. There is more then one possible 
choice for the function x(u). Different choices lead to different possible Lax matrices whose 
eigenvalues lead to distinct linear combinations of the iV conserved quantities. A particularly 
convenient choice is, 

x(u) = — csCKu/in) (10) 

where cs(u) is defined as the ratio of the Jacobian elliptic functions cs(u) = cn(u)/sn(u), 
and K and K' are complete elliptic integrals. The square of x(u) is periodic in u with periods 

2co>i = 2-R-i 

2vrK' , . 

2uj 2 = — = 2niT (11) 

K 
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so 

K' 

K = -*r (12) 
The Lax matrix in the Higgs vacuum equilibrium position has the form 

Lab = (1 - 8 ab )—cs (KUab/iTf) (13) 

7T 

where 

u ab = ^{a-b) (14) 

Conveniently, the matrix L is a circulant, and its eigenvalues are computed by a standard 
formula as, 

N-l 
m=l 

A m = -—cs (2Km/N) 

7T 

z a = e 2ma/N (15) 

Using the peridic properties of cs(u), the eigenvalues can be rewritten in the form of a 
discrete Fourier transform 

2nam 

X a = i 2^ A m sin (16) 

m=l ^ V 

In the continous large N limit this becomes a continuous Fourier transform, and the eigen- 
values are the Fourier coefficients of the function cs(u). Here we note that the Lax matrix 
is manifestly traceless. 



2.1 Semiclassical limit 

Before giving a general formula for the eigenvalues in the Higgs vacuum, we examine the 
semiclassical limit of the integrable system, g — > limit r — > ioo. In this limit potential 
becomes trigonometric, 

v = £ U^ix a -x b ) + 12) (17) 

The equilibrium configuration for the system coincides with the configuration corresponding 
to the Higgs vacuum in the general case. The Lax matrix for this system in the equilibrium 
position reduces to, 

1 ir(a-b) 



J ab — TTl 1 Uab 



1 - 5 nh ) COt 



N 
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with the Lax matrix function x{u) 



1 u 

x{u) = -coth- (19) 

We have again chosen the Lax matrix function x(u) for which L is a circulant, and so is 
convenient for the computation for the eigenvalues. Comparing the form of the Lax matrix 
to the general case (0), the general CM system is recognized as the elliptic deformation of 
the integrable system with 1/ sinh 2 potential (Calogero-Sutherland system). 
The eigenvalues are computed as 

i N-l 

A a = - £ A m z™ (20) 

Z m=l 



where 



7cm 

A m = — cot — — 

2iria/N 



e 2ma/N (21) 



A a can also be rewritten as 



i ^L~l Tim 2nam i , 
X « = "2 ^ cot W sin nv _ = ~2 ( a) a = 1 --- N ~ 1 

m=l 

X N = (22) 

The summation identity used here is the first one in a series of diophantine identitities first 
given in ||. In general these relations are highly non-trivial and conderations of special 
properties of the Lax matrix of the Calogero-Sutherland system actually provided the proof. 
In this particular case the result is not hard to prove by elementary means. In the following 
we will derive elliptic deformations of these identities. 

The above formula for the eigenvalues of the Lax matrix, nearly matches the classical field 
theory formula A^ = — i(N— 2k— 1)/2. In semiclassical field theory, the N eigenvalues of $ are 
evenly distributed along the imaginary axis with unit spacing centered around the origin. In 
(p2|), only N — 1 eigenvalues are distributed in this way and there is an extra zero eigenvalue. 
This discrepancy is suppressed by a power of 1/N and the two eigenvalue distributions agree 
in the large- iV limit. This supports our identification of the eigenvalues of the Lax matrix 
with the field theory VEVs. As explained in the introduction, this identification can now 
be used to yield predictions for the D3 brane distributions in the other vacua which are 
unambiguous in the large- N limit. 

2 An easy way to prove the relation is considering the difference of two consecutive eigenvalues. 
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2.2 General formula for the eigenvalues 

A general formula for the eigenvalues arises using the Fourier expansion of the function cs(u), 



cs(u) = — 



ttu . y-^ q . nixu 

cot 4 > — sin 

2K f-; 1 



2n 

(23) 



n=l 



+ q 2n K 



where 

q = e~ n ^ = e inT (24) 

is the factor associated with a single Yang-Mills instanton. 
Substituting the expansion 



% nm . 2iram ^ q 2n _ 2imm 



A a = — > cot — — sin — — h 2 > > — 

z m=l iV iV m =l„=l -"-Tf/ 



sin e ^^/iv / 25 n 



The first piece in the summation gives back the semiclassical result for the eigenvalues 
(J22|). The remaining terms are instanton corrections. 

Next we evaluate the eigenvalues for the general case keeping the full series in the summa- 
tion. Since the sum is absolutely convergent, the order of the summation is interchangable. 
The double summation then can be performed by rewriting 

oo JV-1 _2n 

K = -UN - 2a) - i £ £ « (e^(^) - e^M) (26) 

z n=l m=l i + ? 

Summing the resulting exponentials, the first summation selects 

n = IN -a I > 1 

n = a -IN I < (27) 



and the sum becomes 



. oo 2(lN-a) 1=0 J2(a-IN) 

1 1 at on ■ at I \ ^ y \ - q 



\ a = — (N-2a)-iN\y - ^ , - V - , , v 

2 V 7 ^ 1 + g2(W-a) ! + g 2(a-UV) 

^ oo / (2l-l)N qN-2a y(2l-l)N „2< 

= --(jV - 20) - ZJV^ I — 9 ( 2J _ 1)Jv ^_ 2a - x - q ( 2 l-l)N q 2a-N J 



The summation is recognized as the expansion of the derivative of logarithmic Jacobi theta 
function 
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The final result of the summation can be written as 



*. = -!(*- 2.) -itfSMg ,30) 



with 



q = q N = e Ni * T 

z a = l -(N-2a)vT (31) 

This is a general formula for the eigenvalues as a function of iV and the effective coupling in 
the Higgs vacuum f = Nr. 



3 Confining vacuum 

In this section we compute the eigenvalues in the unique vacuum of the theory which exhibits 
purely magnetic confinement. 

As reviewed in Section 1, the confining vacuum in question is the image of the Higgs 
vacuum considered in the previous section under the electromagnetic duality transforma- 
tion t — > — i. As the moduli Uk each have holomorphic modular weight k, the individual 
eigenvalues should transform with modular weight one. 

Under the S transformation the Jacobi theta function transforms according to Jacobi's 
formula 

(32) 



I ,2 



Z /-KIT / 







I) 









Os(z\t) = (-it) 2e 

Performing the S transformation including the modular weight we get for the eigenvalues in 
the confining vacuum 

(33) 

where r = is the effective coupling in the confining vacuum, in terms of which q and z a 
are defined as 

q = e~™l f 

Notice that this result depends on r only via the effective coupling constant f = t/N. In 
the large- N 't Hooft limit q = e _7rA . In order to analyze the small A limit, we need to use 
the modular properties of our result in the effective coupling f. In particular, performing 
the transformation f — > — 4 and we find 

r 

~ 0' 3 (-rZa\f) 2tfz a 

e,{z a \-\)~ T e 3 (-fz a \f) + vr 
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Then the eigenvalues can be rewritten as 





~TZ a 


r) 


03 ( 


-fz a 


r) 



(36) 



Here the anomalous piece of the 9 functions under the S transformation exactly cancels the 
leading constant term in A a . In the small A limit f — > zoo, e mr — > 0, and in this limit only 
the first term in the expansion of 9 functions survives. This corresponds to the decoupling 
limit where the theory flows to the 't Hooft limit of the minimal M = 2 SUSY Yang-Mills 
theory with gauge group SU(N). The limiting formula for the eigenvalues is 



A a = 2sin(2f5 a )e i 



-2e _7r / A si 



sm 



2air 

~n~ 



(37) 



In this limit, the eigenvalues computed from the Lax matrix can be compared with the 
results of Douglas and Shenker Jl5| for the minimal Af = 2 theory. In the large- iV limit we 
can find the density p of the eigenvalues normalized to lie between zero and one. The result 
is the Wigner distribution 



P(A) 



1 



X'(x) 



x=x(\) 



7T 



(38) 



This agrees with the large-iV limit of the distribution found by Douglas and Shenker as 
computed in JTBI PI. 



4 General vacua in the supergravity limit 

To compute the eigenvalue distributions in the supergravity limit, we first examine the limit 
in the confining vacuum. The eigenvalues in the confining vacua are given in (^), and in 
terms of iV and the 't Hooft coupling A can be explicitly written as 



A a = A 



1 a 00 

---+y 

2 N 



e -39TA(i-#) 



e -27rA(/-l+f ) 



1 _|_ e -27rA(Z-f ) i + e -2n\(l-l+§) 



(39) 



The supergravity limit is taken iV — > oo, while keeping the 't Hooft coupling A fixed, 
large. Taking the N — > oo limit ^ becomes a continuous parameter ranging from to 1. In 
this limit we have the continuos distribution 



X(x) = A 



2"* + E 
z 1=1 



! Our normalization is such that J Q p(x)dx 



1 _|_ e -2nX(l-x) 
= 1. 



I _|_ e -27rA(Z-l+x) 



(40) 
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When A is large the summation terms are exponentially supressed, so in the supergravity 
limit the the confining vacua are 



X(x) =AQ-z) (41) 

corresponding to a constant distribution. (We note that at the endpoints x = and x — 1 
there is a constant shift in the eigenvalues. A similar endpoint effect can be seen when ana- 
lyzing the eigenvalue distribution from the Seiberg-Witten curve. Such a constant endpoint 
shift does not effect the constant distribution of the eigenvalues.) 

It is convenient to normalize the uniform distribution to range one, so in the following 
we will consider the suitably normalized eigenvalues. In the confining vacuum this means 
dividing by the 't Hooft coupling A, so the distribution we arrive at the supergravity limit is 

A» = \ - x (42) 

where the normalized eigenvalues are denoted by A. 

Next we examine the supergravity limit of the eigenvalues in a general vacuum. The set 
of vacua are labelled with three integers (p, q, I) with the coupling constant 

f = pq = N l = 0...q-l (43) 

Q 

The Higgs vacuum corresponds to (N, 1, 0) while the confining vacuum is (1, N, /). An 
SL(2, Z) transformation performed on a generic vacuum takes it into another vacuum, and 
also changes the coupling constant to a different coupling f. Then there exists another SL(2, 
Z) transformation which maps this vacuum to one element of the sets of the original vacua, 
characterized with different integers (p', q', I'). In this way SL(2, Z) transformations map the 
set of vacua among themselves. 

A general vacuum is th then connected to the confining vacuum by a modular transfor- 
mation with weight 1: 

*M = -q-;JW f 2 ^) < 44 > 

cr + a \ct + a J 

The weight can be seen for example from the description of vacua from the curve, in which 
the Riemann zeta functions transform with modular weight 1. 

In a general vacuum the supergravity limit is taken by taking N to infinity, while keeping 
the dual 't Hooft coupling fixed and large. For a general vaccum connected to the Higgs 
vacuum by a modular transformation (03) we find in the large N fixed Xd limit 



I oo / e -27rA D (i-x) e -2Tr\ D {l-l+x) \ 

K X ) = 2~ X + + e -2n\ D (l-x) ~ I + e -2w\ D (l-l+x) ) ( 45 ) 
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where Xd is defined from td 



ar+b 
cr+d 



as Ad = IN/tb- Here we normalized the eigenvalues 



appropriately by Xd/(ct + d). The modular factor comes as the weight of the weight-1 mod- 
ular transformation. In the supergravity limit the summation is exponentially suppressed, 
and the eigenvalues for the general vacua are 



1 

2~ X 



(46) 



In conclusion, in the supergravity limit we find a uniform distribution for a general 
vacuum. The corrections to the supergravity limit are the summation terms with the ex- 
ponential suppression in (^). These terms represent instanton contributions to the uni- 
form supergravity distribution coming from (m, n)-strings which wrap the two-sphere in the 
Polchinski-Strassler geometry. The contribution of instantons can be rewritten by expanding 
the denominators and collecting the powers of each exponential term e ~ m . In this way 
we arrive at the alternative expression 



1 



-*+E 



m=l 



( 



E 



E 



d\m 
-k mod JV 



d\m 
d=k mod JV 



/ 



(47) 



This form of (|4"5| ) exhibits the coefficients the instanton sum explicitly. 



5 Condensates 

The physical quantities which are computed from the eigenvalue distribution are the conden- 
sates Uk = Tr$ fe . These essentially correspond to the integrals of motion for the integrable 
system, which are expressed as the symmetric polynomials of the eigenvalues of the Lax 
matrix. 

In this section we compare some of the condensates computed from the eigenvalues to 
earlier results for these. First we examine the semiclassical limit of the condensate u<i = Tr$ 2 
in the Higgs vaccum. 

The precise connection between U2 and the symmetric polynomial constructed from the 
squares of the eigenvalues is encoded in the Lax equations. In order the Lax equations to 
give Hamilton's equations, the square of the function x(u) in the Lax matrix is related to 
the potential V(u) by a constant, as described in In the semiclassical limit 

Vab = V{u ab ) = 1 + — 

4 sinn ±Uab lz 
x ab = X(u ab ) = - coth \u ab (48) 



so in this limit 



2 2 



C=-\ (49) 
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We also have 

h= 1 -TrL 2 = J2xl b = E(^-C) (50) 

a>b a>b 

from which 

u 2 = Y,v ab = h 12 = 2 ^ a 12— (51) 

a>b a=l 

Substituting the semiclassical result for the eigenvalues (p2| ) 

N(N 2 - 1) 

U2 = (52) 



in perfect agreement with the semiclassical result of jl3 . 

Next we compute the condensate u 2 in the confining vacuum. From the eigenvalues in 
the confining limit (^) we find 

i N N 
« 2 = h = ? E A ' = 2e" 2 x 2 sin 2 — - = Ne~- = Ne~ (53) 

/ o=l 0=1 iV 

This is in agreement with the result for u 2 computed in [|l3j in terms of Eisenstein series 

U 2 = ^(e 2 {t)-^E 2 (?-\) (54) 



24 V y ' N \N 
upto a vacuum independent shift function. The relevant part of U 2 is the second piece 



U' 2 = -—E 2 ( — ) ^ Ne~ (55) 



24 ViV 

in agreement with u 2 computed from the eigenvalues. 
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A Distribution of vacua from the Seiberg-Witten curve 



The Seiberg-Witten curve of the Af = 2* SU(N) theory can be obtained from IIA/M-theory 
brane construction |TU| where one considers a brane system suitably constructed from a 



single NS5 and N D4 branes. The Coulomb branch is obtained as a branched N-fold cover 
of the torus with complex structure r 

4zvr 9 iN 6 
T = ^ + ^ = - + 2^ (56) 

where r is the complex gauge coupling, and A = g 2 N/4ii is the 't Hooft coupling. 

The massive M = 1* vacua are the singular points on the Coulomb branch where the 
Seiberg-Witten curve maximally degenerates. At these points, the curve becomes an un- 
branched cover, a torus itself, with complex parameter 

f = pr±j_ pq = N i = Qq _ 1 (57) 
Q 

A generic vacuum is thus characterized by the 3 integers (p,q,l). The Higgs vacuum cor- 
responds to (N, 1, 0) while the confining vacuum is (1, N, I). For the confining vacua, it 
is enough to consider the case I = 0, because I nonzero is connected to this by a modular 
transformation. 

The general SW curve can be described in terms of coordinates of the D4-branes as 

II(« -«.(*)) = (58) 

a=l 

here z is on the covered torus with complex structure r, which has periods r = io%ju)\ . 

The M = 1* vacua correspond to the singular points where the curve maximally degen- 
erates to a torus with complex structure f. The coordinates at these vacua, taking the case 
1 = 0, are given as E 



-m 



v sr ( z ) = NmQ{z — Z\ + 2sui + 2ru> 2 ) 

fp-lq-l _ 

J2J2(( Z + 2tuJ i + 2uuJ 2) + 2ps((pui) + 2qr((quj 2 ) I (59) 



t=0 u=0 



where s = 0, l..q — 1 and r = 0, l..p — 1. ( is the Weierstrass zeta function with complex 
structure f and periods LO\ = quii, and uj 2 = pou 2 , where we fix ui\ = in. Z\ denotes the 
position of the NS5 brane, and m is the hypermultiplet mass. 

In the following we determine the distribution of Higgs and confining vacua in the large 
N limit. 

In the Higgs vacuum (N, 1, 0) we denote the r = 0, 1..N — 1 coordinates by vf(z) 



?(z) = Nm((z - z x + 2rinr) -m(^({z + 2umT) + 2r((NinT)\ (60) 

\u=0 / 
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We are interested in the distribution of vacua around the their average, so we redefine the 
coordinates to be centered 

V r h (z)=v?(z)-v H (z) (61) 

so that 

N-l 

E V r B {z) = (62) 

r=0 

The new coordinates are then expressed as 



N-l 



V r h (z) = Nm((z -zx + 2rmr) - m ^ ((z - z l + 2umr) + m(N - 1 - 2r)C(Ninr) (63) 

Taking the limit N large while keeping the 't Hooft coupling A fixed, r — > ioo. The ( 
function is periodic with periods 

cui = quj\ = in 

U>2 = P^2 = NlTlT (64) 



and its complex structure is 

f = 



t = — = Nt — > too (65) 



In this limit the ( function reduces to 



12 2 2 

and 

T r h , s Nm , ( z — Zi . \ 

V r k (z) = — COth \—^- L +TI-KT J 



C(z) = -— + icoth- (66) 



m ^-^ i ( z — z i ■ \ m ,*r -, _ , Ninr 

coth I + uvkt I + — (N - 1 - 2r) coth (67) 

2 u=0 V 2 / 2 2 

In the r — > ioo limit this becomes 

V r h (z) = mr — coth r > 

Zi Zi z 

V h (z) = mr + m(JV " 1} coth^l r = (68) 

z z 

The fact that the r = coordinate is different has no significance in the continuous iV — > oo 
limit, since this is of measure zero. 
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To determine a distribution for the coordinates, we scale them between and 1, that is 
divide by N. Defining the parameter 

x = ^ (69) 

and 

V h (z) 

Kx) = (70) 

we find 

h(x)=mx — m/2 (71) 

Here the z-dependent piece is which is of order l/JVis dropped. Rescaling the coordinates by 
m we find the constant distribution in the Higgs vacua in the large N fixed 't Hooft coupling 
limit 

Pu{h) 



h'(x) 



. = 1 (72) 

x=x(h) 



Here we denoted the rescaled coordinates by h. (The mass m in the integrable system 
description is absorbed in the scaled variables in which the superpotential is written. ) 

The constant distribution computed from the curve agrees with the constant distribution 
found in terms of the eigenvalues of the Lax matrix in the semiclassical limit fl22[). 

One can similarly analyze the supergravity limit of the confining (1, N, 0) vacua, with 
the centered coordinates 



N-l 



t=0 



V s c (z) = Nm((z - zt + 2sm) - m ^ CO - *i + 2ti?r) + m(N - 1 - 2s)((Nm) (73) 
C now has the periods 



uji = Nm 

uj 2 = «7rr (74) 



and complex structure 



T = — = ~ T 75 
oj\ N A 

Here the 9 term in the gauge coupling is neglected. In the large TV limit with fixed large 
't Hooft coupling A (the supergravity limit) f — > 0. In order to get an expansion of the 
( function, a modular S-transformation has to be performed. The modular transformation 
acts on ( as 

. ar + b 

r = ; 

cr + a 

C(z\t>) = (cr + d)((z(cT + d)\r) (76) 
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Performing the S transformation on ( when A is large, so in the supergravity limit, the £ 
function can be replaced by its semiclassical limit. Then the supergravity (large N, large 
A) limit can be similarly analyzed as for the case of the Higgs vacuum. We find a constant 
distribution in agreement with the distribution computed from the integrable system fl46|). 
Because of the modular properties of the ( function, the coordinates of the Higgs and con- 
fining vacua as well as their large N distributions are related by the modular transformation 
in r = iN/ A as 




(77) 
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